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Abstract 



We study series of the stationary solutions with asymptotic flatness properties in the 
Einstein-Maxwell-free scalar system because they are locally equivalent with the exterior 
, solutions in some class of the scalar-tensor theories of gravity. First, we classify spherical 

^ ■ exterior solutions into two types of the solutions, an apparently black hole type solution 

' and an apparently worm hole type solution. The solutions contain three parameters, 

and we clarify their physical significance. Second, we reduce the field equations for the 
' axisymmetric exterior solutions. We find that the reduced equations are partially the 

same as the Ernst equations. As simple examples, we derive new series of the static, 
axisymmetric exterior solutions, which correspond to Voorhees's solutions. We then show a 
non-trivial relation between the spherical exterior solutions and our new solutions. Finally, 
since null geodesies have conformally invariant properties, we study the local geometry 
of the exterior solutions by using the optical scalar equations and find some anomalous 
behaviors of the null geodesies. 



1 Introduction 



Recently, as natural alternatives to general relativity, the scalar-tensor theories have been 
studied by many theoretical physicists. In these theories, the gravity is mediated not only by 
a tensor field but also by a scalar field. Also, such theories have been of interest as effective 
theories of the string theory at low energy scales . 

Several theoretical predictions in the scalar-tensor theories have been obtained (see e.g. 
Ref.§~§]). It has been found that a wide class of the scalar-tensor theories can pass all the 
experimental tests in weak gravitational fields. However, it has also been found that the scalar- 
tensor theories show different aspects of the gravity in the strong gravitational fields in contrast 
to general relativity (see e.g. Ref. (5|). It has been shown numerically that nonperturbative 
effects in the scalar-tensor theories increase the maximum mass of an isolated system such 
as a neutron star [2-5]. In these works, numerical methods play an important role, and, in 
interpreting the numerical results, a static, spherically symmetric vacuum solution (a spherical 
exterior solution) is matched to the numerical solution. 

The spherical exterior solution is an analytic exact solution, and, as far as the authors 
know, only few exact solutions have been known. In particular, the axisymmetric exterior 
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solution corresponding to the Kerr solution must be of significant interest. Motivated by this 
fact, we study series of the stationary solutions with asymptotic flatness properties in the 
Einstein-Maxwell-free scalar system because, as summarized in Appendix A, they are locally 
equivalent with the exterior solutions in some class of the scalar-tensor theories of gravity. The 
field equations are given by 



where T^ v is a energy- momentum tensor of the Maxwell field (see also Appendix A). 

In this paper we first classify spherical exterior solutions into two types of the solutions, 
an apparently black hole type (ABH) solution and an apparently worm hole type (AWH) 
solution. We then clarify physical significance of the parameters contained in the solutions 
and show that the ABH solution and the AWH solution correspond to, respectively, the mass 
dominant case and the charge dominant case of the Reissner-Nordstrom solution. Second, we 
reduce the field equations for the axisymmetric exterior solutions and show that the reduced 
equations are partially the same as the Ernst equations. As simple examples, we derive new 
series of static, axisymmetric exterior solutions, which will be referred to as scalar-tensor- Weyl 
solutions. We then show a non-trivial relation between the spherical exterior solutions and 
our new solutions. Finally, we study the local geometry of the exterior solutions by using the 
optical scalar equations and find some anomalous behaviors of the null geodesies. 



2 Spherical exterior solutions 

2.1 Classification of the exterior solutions 

In this section, we derive static, spherically symmetric exterior solutions with the electric 
field, which are hereafter referred to as spherical exterior solutions. In the Just coordinate, the 
metric becomes 101: 



where 7 and A are functions of x- Non-vanishing components of the electromagnetic tensor, 
F^, are Fqi = — F10 = E(x)- The Maxwell equation is reduced to 



Rf+v = SnTf,,,, + 2cp l/t ip l 
□^ = 0, 



(1.1) 



ds 2 = -e^dt 2 + e-~<d X 2 + e^dti 2 , 



(2.1) 




(2.2) 



(2.3) 



^( 7 " + A' 7 V - Q 2 e 2 ^ A \ 

X - (7" - 2A" - A' 2 - 7 ' 2 + A'7 V - 2</ V - Q 2 e 2 ^~^ , 

e^ A + i(7"-A" + AV-A ,2 )e^ = Q 2 e 2 ^" A \ 

(ip" + AW = 0. 



(2.4c) 
(2.4d) 



(2.4b) 



(2.4a) 
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With ( pT4a| ) and ( fOc] ), one obtains 



(e A )" = 2. 



(2.5) 



One finds two types of the solutions to ( |2.5| ), which are referred to as an apparently black hole 
type (ABH) solution and an apparently worm hole type (AWH) solution, respectively: 



ABH 
AWH 



X -ax 
9 a 2 

x 2 + ^ 



( X >a>0), 
(-00 < \ < oo), 



where a is an integration constant. 

When a^O, the equation (2.4d) is integrated as 



ABH ip = ipo + 
AWH if = ip + 

where fo an d c are integration constants [2-5]. 
The equation (2.4a) is integrated as 



c . 


I a\ 


-hi-- 


a 


V X) 


2c 


\ c> 




arctan - 


a 





(Ye 



4Q 2 e 7 + e6 2 , 



(2.6a) 
(2.6b) 



(2.7a) 
(2.7b) 

(2.8) 



where b is an in te grati on con stant , and e is 1 in the ABH solution and —1 in the AWH solution. 
With fl2.4a| ), Q2.4b|) and ( $A<\) , we obtain 



A — 7 



-A 



V —Qe 



2 o7 -2A 



which is reduced to the following relation among the integration constants: 

ABH a 2 - b 2 = Ac 2 — ► a > b > 0, 

AWH b 2 - a 2 = Ac 2 — > b > a > 0. 



Since e 7 — ► 1 as \ — > 00 m the asymptotically flat spacetime, with (2. 
following inequality for the AWH solution: 

b 2 < 4Q 2 . 



(2.9) 

(2.10a) 
(2.10b) 

we obtain the 
(2.11) 



That is, the AWH solution is charge dominant, and no AWH solution exists when Q = 0. 



2.2 The ABH solution 

The metric function, 7, in the ABH solution is obtained as follows: 

(2.12a) 



B ( b x 2 



(1-5) 
where 



B^k[X- a -\\ k^^/^ll^-' <1. (2.12b) 

V4Q 2 + b 2 + b 



3 



In the limit, Q — > 0, the solution coincides with the previously known solution [2-5]. In the 
limit, b — ► 0, the solution is reduced to the following form: 



e' = 



1 — — In I 1 — — 

X 



(2.13) 



When a = b, i.e., c = 0, the Just coordinate, x> is related to the Schwarzschild coordinate, 
r, as 



The metric function, 7, becomes 



kb 
1-k 



1 - k 



X = r 



kb 
1-k' 



(2.14) 



where 



1-k' 



1-— ) (I-—) =1- — + (2.15a) 



(2.15b) 



kb 1 + k 
-, 2m = - 

1-k 1-k 



b=^4Q* + b 2 . 



With ( 2.15b ), we obtain the following inequality: 

< b 2 = 4{m 2 - Q 2 ). 



(2.16) 



That is, the ABH solution corresponds to the mass dominant (m 2 > Q 2 ) Reissner- Nordstrom 
spacetime and has two event horizons, r + and r_, when a = b. 

When a > b, the Just coordinate, x, is related to the Schwarzschild coordinate, r, as 



a\(l-B) 2 (Q 
X 1 

X 



B 



(2.17) 



which vanishes at x — a i an d Xb = o,/ (l — fc _a / b ) < 0. Moreover, we find that x — a : X = 
and x — XB, correspond to, respectively, r = r + , r = r_ and r = when a — b and that the 
null surface, x — a i becomes a singularity when a > b. 



2.3 The AWH solution 



When a > 0, the metric function, 7, in the AWH solution is obtained as 

2 

4 \Q, 

where the constant, [3, is defined by 



Cretan 

a \ a 



■0 



nb 



4Q^ 
b 2 ' 



(2.18a) 



(2.18b) 



As mentioned previously, the AWH solution does not exist when Q = 0. 
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When a = 0, the AWH solution becomes 

e A = X 2 , 
<P = <A> - 



X 



of = 



sec 2 P 



where 



sec (3 — 



Q 2 



When a = b, the AWH solution is reduced to the following: 



, i 2 * 
arctan — — 





2 b 2 


e A = 




e 7 = 










<P = 





where 



sec' 



4g^ 

6 2 ' 



Note that (2.22) is equivalent with the following: 



b / b 2 

sin/3 = ± — — , cosp = < / 1 — — — j 



Then (2.21b) becomes 



4Q 2 



"2Q 



1 + 4— 

6 2 



4g 2 



L _ X 

4g 2 + q, 



and the Just coordinate, \i is related to the Schwarzschild coordinate, r, as 



r = e 



(A-7)/a _ 



fo 2 



When we take r , we have 



6 2 1 O 2 

-l + 2 x lQ 2 - U -- + \ = l- — + \, m. 

4 7* 7* 7* 7* 



2m Q 2 



Q 2 - ~ < 0, 



(2.19a) 
(2.19b) 

(2.19c) 
(2.19d) 

(2.20) 

(2.21a) 

(2.21b) 
(2.21c) 

(2.22) 
(2.23) 



(2.24) 



(2.25) 



(2.26) 



which is a negative mass solution and should be discarded as an unphysical exterior solution. 
When we take r + , we have 



1 - 2 



b 2 1 , Q 2 _ i 2?7i ( Q 



... 2 / 6 2 

-— - + ^- = 1 + ^-, 0<m=WQ 2 -— . 



4 r 



(2.27) 
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That is, the AWH solution corresponds to the charge dominant (m 2 < Q 2 ) Reissner- Nordstrom 
spacctime and has a naked timelike singularity, r — 0, when a = b. 

When a < b, the Just coordinate, Xi is related to the Schwarzschild coordinate, r, as 



X + -j | cos 



-arctanl^ 



(2.28) 



Note that r = at x — Xw, where 



Xw = ^ tan 



6 V2 







(2.29) 



The AWH solution has a naked timelike singularity at x = Xw- That is, the AWH solution 
does not have a worm hole, and the asymptotic region defined by x ~~ * — 00 does not exist. 



3 The Ernst equations in scalar-tensor theories 



In the stationary, axisymmetric spacctime, there exist two Killing vectors, £ = dt and 
T] = d<p. It is shown that the energy- momentum tensor, T a @ , of the Maxwell field satisfies the 
following relations Q: 

CT a l(3 Cv S] = V a T a l(3 Cr] S] = 0. (3.1) 

Moreover, one finds that the similar relations hold for t}$ = 2d^ipd v ip. Accordingly, with the 
theorem 7.1.1 in Ref.]7j], the metric can be reduced to the following form: 

ds 2 = -e 2 ^(dt - tod<P) 2 + e' 2 ^ [e 2 ^(dp 2 + dz 2 ) + p 2 d<p 2 } , (3.2) 

where the metric functions, ip,^y and u>, are functions of x 1 = p and x 2 = z. 

After long and complicated calculations, we find that the field equations are reduced to the 
following: 



e 2 ^ 



i (3.3a) 

e 2 ^V 2 £ = (V£) • [(y£) + 2$V$] , (3.3b) 

e 2 ^V 2 $ = (V$) • [ye) + 2l>V$] , (3.3c) 

7.1 = p [(^ a ) 2 - (^ 2 ) 2 ] - 1 [KO 2 ~ K2) 2 ] e 4 ^ 

-p($ i$ i - $ 2 $ 2 )e~ 2V ' + p [(^,i) 2 - (^, 2 ) 2 ] , (3.3d) 

7.2 = 2 pip, lip, 2 - -^-oj.xuu^ - p($, i$2 + "J^^i)^ 2 ^ + 2(xpitp t2 , (3.3e) 

2p 

VV = 0, (3.3f) 
where, for any functions, / and h, 

72f _i3 f_df\ , a 2 / „ u _dfdh t dfdh t _ df 



^ p dp \ dp ) dz 2 1 ^ dp dp dz dz' dx % ^ ^ 



G 



In Ref. ||, one will find the explicit forms of the metric function, u, and the Maxwell field in 
terms of the complex potentials, £ and <&. Note that the first, second and third equations 
are the same as those in general relativity and are referred to as the Ernst equations of the 
Einstein-Maxwell system. The remaining equations contain the scalar field contributions, and, 
as for the metric functions, effects of the scalar field only appear in 7. 



4 Scalar-tensor- Weyl solutions. 
4.1 The reduced Ernst equation 



In this section, we shall consid er st ati c, ax isymmetric vacuum solutions. In this case, 
oj = $ = 0, and the field equations ( 3.3a )~( 3.3f ) are reduced to the following: 

3 2V 



e"" r = B, 
72 < 



£V Z £ = (V£) • (V£), 

7.1 = P [(ti) 2 -(t2f]+p[M 2 -M 2 ] 

7.2 = 2pV,ll/\2 + 2/5(^1(^2, 

vV = 0, 



where £ is a real function. One finds that the Ernst equation (4.1b) is reduced to 

VV = 0. 

That is, <p and ip are harmonic functions. 

We introduce oblate and prolate coordinates, (x,y), defined by 



p = <7^{x 2 + e){l-y 2 ), 
z = crxy, 



(4.1a) 
(4.1b) 
(4.1c) 
(4.1d) 
(4.1e) 

(4.2) 



(4.3a) 
(4.3b) 



where a is a positive constant, and e = ±1. The cases, e = 1 and e = — 1, are referred to as 
oblate and prolate, respectively. The Ernst equation is then reduced to 



d_ 

dx 



(x 2 + e) 



dip 
dx 



d_ 

dy 



2^ 



dy 



0. 



(4.4) 



The similar equation holds for ip. 



4.2 Prolate solutions 



The simplest prolate solution for -0 is given by 



1> = iu(?-±) 1 (4.5) 



2 \x + 1 

where 5 is an integration constant. The similar solution for tp is given by 



d fx — 1\ 
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where (po and d are integration constants. Then the corresponding solution for 7 becomes 

A 2 



.27 



X Z -1 

x 2 — y 2 



where 



A 



(4.7) 



(4.8) 



For completeness, an explicit form of the metric is shown below: 



ds 2 



x-l 
x+1 

x 2 -\ 
x 2 — y 2 



dt 2 + a 2 



x-l 



x + 1 



—8 



dx 2 



dy 2 



x 2 — 1 1 — y 2 



+ {x 2 -l){l-y 2 W 



(4.9) 



Note that, though ip and ip are y-independent, 7 depends on both x and y. The metric 
contains two parameters, 8 and A. Now we show that the metric is reduced to the previously 
known one when we take specific values of the parameters. When A = S, i.e., d = 0, the metric 
is reduced to one of the Weyl series of the solutions in general relativity (Voorhees's prolate 
solution ||). Therefore, we refer to (4.E) as the scalar-tensor- Weyl solution. Voorhees's prolate 
solution contains the Schwarzschild solution as a specific case, S = 1, and the Schwarzschild 
coordinates are related to {x, y) as 



x = 1, y = cos9, a = m, 

m 



(4.10) 



where m is a usual mass parameter. 

The case, A = 1, is the most interesting one. We introduce new parameters, a, b and c, 
and new coordinates (the Just coordinates) by 



and 



, 2c a b 

d = — , a = -, d = -, 
a 2 a 



x = 1 , y = cos ( 



(4.11) 



(4.12) 



Then we find that the solution coincides with the spherical exterior solution with a vanishing 
electric charge: 



ds 2 



dt 2 



d X 2 



x 2 li 



sin 



(4.13) 



Moreover, when 5 = 1, i.e., a = b, the solution is reduced to Schwarzschild's one. 

The right hand sides in ( 4.1cj ) and (4. Id) are in the forms of the superposition of the 
contributions of ip and ip. It is important to note that the superposition of the effectively non- 
spherical i/j-terms corresponding to Voorhees's prolate solution with 5^1 and the spherical 
configuration of tp reduces to the solution containing the spherical exterior solution. When the 
spherical scalar field configuration is added to the -0-terms corresponding to Voorhees's prolate 
solution with 8 = 1, namely the Schwarzschild solution, we cannot have spherical solutions 
unless a = b. The relation among these solutions is schematically shown in Figure 1. The 
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Figure 1 : A schematic picture of series of the scalar-tensor- Weyl solutions. 



scalar-tensor- Weyl solution generically has a singularity, x = 1, whose topology is S 2 , and 
qualitative features of the singularity may be similar to those in Voorhees's prolate solution. 
When A = 1, the singularity, x = 1, becomes a point, which is given x = a in the Just 
coordinate. Moreover, when A = S = 1, the topology of \ = a, corresponding to the event 
horizon, r = 2m, becomes S 2 . 



4.3 Oblate solutions 



The simplest oblate solution for ip is given by 

ip — S ■ arccotx, (4-14) 
where S is an integration constant. The similar solution for <p is given by 

<P = <Po + d ■ arccota;, (4-15) 
where ip and d are integration constants. Then the corresponding solution for 7 becomes 



■ 2 2 \ A 2 

,2 7 _( x 2 + V 2 



e ' = 



x 2 + 1 



(4.16) 



where A 2 = 5 2 + d 2 . For completeness, an explicit form of the metric is shown below: 

ds 2 = _ e 2<5 - arccote ^ 2 + cr 2 e -25.arccot 2; x 

I- A 2 

dx 2 dy 2 



x z + 
x 2 + l 



1 i-r 



{x 2 + l){l-y 2 )d<t> 2 



(4.17) 



When A = 6, the metric is reduced to Voorhees's oblate solution. In contrast to the prolate 
case, the oblate solution is not reduced to the spherical exterior solutions for any S and A. 
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5 Local behaviors of null geodesies 



5.1 The spherical exterior solution 



In this section, we examine local behaviors of irrotational null geodesies in the spherical 
exterior solution with a vanishing electric charge. The metric is rewritten in the Schwarzschild 
coordinate as 

ds 2 = (l -^f dt 2 + (l ^1 VV +rW, 



X(r)j V r ^ (5 i) 

where the mass function, m(r), can be interpreted as a local energy []l0|| . One immediately 
finds that m(r) becomes negative when (a + b) 2 /4b > x > a - Therefore, one may expect that 
there will be significant difference not only in the global causal structures but also in the local 
geometry, depending on whether a is equal to b or not. However, we will find that it is not 
the case and that local behaviors of null geodesies significantly depend on whether a is greater 
than 2b or not. 

Due to the conformally invariant nature of null geodesies, we consider null geodesies in 
the spacetime with a metric, g = A _2 g (see Appendix A). Let k" — dx^/dv be null geodesic 
associated with g with the affine parameter, v. Hereafter, a hat denotes geometrical quantities 
associated with g, and the corresponding geometrical quatities associated with g are denoted 
without a hat. Then the null geodesies associated with g are given as 

dr" » - ~ - 

d\ = A- 2 dv, = — = A 2 k", -> k a W a k" = k a \J a k" = 0. (5.2) 
dX 

Since the local geometrical nature of the spacetime is described by the Riemann curvature, it 
is important to examine optical scalars of the irrotational null geodesic congruence defined by 
Sachs jnj. Let {E( a ) : (a = 1 ~ 4)} = {k, rh, t,t} be a null tetrad such that 

k a V a E{ a) = 0, faEfcEfa = -2&\ a 5 2 b) - 25f ^. (5.3) 

Then the optical scalars, namely, the expansion, 9, and the complex shear, cr, are defined as 

= V^kJH", a = VpkjH", (5.4) 

The optical scalar equations are also conformally invariant, and we define the following quan- 
tities associated with g: 

6 = A 2 8-k fi V»lnA, a = A 2 a, f = At" . (5.5) 
Then the optical scalar equations become |TT], Q 



cr| 2 = —R^k" = U, 



+ (? + \a\ 2 = -- 

dX 2 (5.6) 

^ + 20a = -R^pkW*^ = -C^k^k^tP = F, 
dX 



where C^avp is the Weyl curvature. 



10 



First, we examine circular orbits of null geodesies in the Just coordinate. Without loss of 
generality, we consider null geodesies on the 8 = tc/2 plane. Then the geodesic equation is 
reduced to 



(2b-a) 

L 2 ( a 



x 2 = E 2 -^(l--) = E*-V( X ), 



X V X 



(5.7) 



where a dot denotes a derivative with respect to the affine parameter, A, and E and L are 
integration constants. The circular orbits are determined by the conditions that x = and 
dV/dx = 0. We find that the circular orbit is obtained as x — Xc = (a + 2b)/ 2 for a < 2b, and 
that, surprisingly, no circular orbit exists for a > 2b in contrast to the case in the Schwarzschild 
spacetime. 



Dyer [|12| has obtained general solutions to the optical scalar equations (5.6) in the general 
static, spherically symmetric spacetime, and we summarize his results in Appendix B. In the 
static, spherically symmetric spacetime, one finds that the Weyl driving term, F, and the 
shear, a, can be regarded as real quantities without loss of generality p3 |. By defining the new 
optical scalars, C±, as 



dX 

Dyer has obtained the following equations 

d 2 C ± 



d lnC*± =6±a, (5. 



d\ 2 

We evaluate the Ricci and Weyl driving terms, TZ and F, as 



(a 2 - b 2 ) 

TZ = - 



(K±F)C±. (5.9) 




. (,_£)" (*)», (5.10) 



h 



2 



a 



(2b— 3a) 



F= 4x E V~x) {^X-(a + b)(a + 2b)}, 

where h = L/E is an impact parameter, and the affine parameter, A, is chosen such that E = 1. 
It is immediately found that the Ricci driving term, TZ, vanishes at the perihelion where k 1 = 0. 
A particularly interesting case is the circular orbit on which k 1 = identically. If F vanishes on 
the circular orbit, the image shape of an infinitesimal light ray congruence remains unchanged. 
The condition, F = 0, is reduced to x = Xf = (a + b)(a + 2b) /(6b) for a > 2b. When a < 2b, F 
is strictly positive. Since the circular orbit exists only when a < 2b, the condition, TZ = F = 0, 
is never satisfied on the circular orbit. However, the condition, TZ = F = 0, can be satisfied at 
X = Xf for the scattering orbit when the impact parameter, h, is chosen as 



{a + b)(a + 2b) 
h= 6b 



(a-b)(a-2b) 
(a + b){a + 2b) 



TZ(xf) = F( X f) = 0. (5.11) 
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That is, the gravity effectively vanishes at \f in this sense. 

In summary, we have found that local behaviors of null geodesies significantly depend on 
whether a is greater than 2b or not. When a < 2b, the behaviors are similar to those in the 
Schwarzschild spacetime. That is, the circular orbit exists, and the Weyl driving term, F, is 
strictly positive. When a > 2b, on the contrary, there is no circular orbit, and F changes its 
sign at xf- 

These behaviors of null geodesies can be compared with those in the Reissner-Nordstrom 
spacetime: 

= - (i - ^ + * 2 + ^ - v + 7?r ^ + r w - (5 - i2) 

The Ricci and Weyl driving terms are calculated as 

Q 2 h 2 3h 2 (m Q 2 \ , s 

1Z = -^-, F= — (--\). 5.13 

r 4 r 4 \ r r z J 

One immediately finds that F vanished at r — rp = Q 2 /m. However, when m 2 > Q 2 , rp is 
inside the event horizon, r + = m + y/ m 2 — Q 2 . When 9m 2 > 8Q 2 , two circular orbits exist: 

r = r± = - (im ± yj%m 2 - 8Q 2 ) . (5.14) 



> Q 2 , one finds that r% > r + > rp > r c > r_, where 



r_ = m — 



y m 2 — Q 2 . A particular case that rp — r c occurs only when m 2 — Q 2 , and one finds that 
rp = r^ = rn = r± and that rf, — 2m. 



General solutions to (5.9) with (5.10) can be obtained by Dyer's formula (Appendix B.l), 
however, we have only examined quantitative behaviors of the Ricci and Weyl driving terms, 
1Z and F. Further quantitative studies of the solutions and their astronomical applications, 
especially the implication to gravitational lens effects, will be discussed in the forthcoming 
paper. 



5.2 The scalar-tensor- Weyl solution 



It is difficult to study generic null geodesies in the scalar-tensor- Weyl solution. Therefore, 
we shall examine specific null geodesies on the y = plane. When y = 0, the geodesic equation 
is reduced to 

1 / X — 1 \ • _ , n / X — 1 



E = - ■ ) t, L = (x 2 - 1) 



a \x + 1 J \ x + 1 

where W is a function of x defined by 



,r 2 



dx 



dW \x 2 -l 



(5.16) 



a dot denotes a derivative with respect to the afline parameter, A, and E and L are integration 
constants. The circular orbits are determined by the conditions that W = and dV/dx = 0. 
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Wc find that these conditions are reduced to x = x c = 28. Since x should be larger than unity, 
the circular orbit exists when 8 > 1/2. 

We evaluate the Ricci and Weyl deriving terms, 7Z and F, (see Appendix B.2) and find that 



F = 



\2 



n = — 



i 



a 2 x 2 (x 2 
A 2 - S 2 



1) v^ 2 - 1 



A 



(x 2 - 1) : 



r(fc 



1\2 



1 



a 2 x 3 (x 2 — 1) \x 2 — 1 



6(A 2 



x 2 - 1 



x 

x - 1 



x - 



1 



(h/af 



1 



2(5' 



2 a 



(5.17) 



{<5(A 2 - 1) - (A 2 + 28 2 )x + Z5x 2 } 



where h = L/E is an impact parameter. One finds that they are reduced to ( 5. 10| ) when 
A 2 = 1. We show K and F in Fi gure 2. We numerically find that the qualitative features of 1Z 
and F are the same as those in the spherical case (A = 1) and that the significant qualitative 
changes in 1Z and F are seen depending on whether 8 < 1/2 or not. 



6 Discussions 



After summarizing the spherical exterior solutions in the scalar-tensor theories of gravity, 
we have derived a new series of the axisymmetric exterior solutions, which are reduced to 
Voorhees's solutions in general relativity. We have found that the prolate solution, namely 
the scalar-tensor- Weyl solution, contains the spherical exterior solution as a special case in 
which we take some non-trivial choice of the parameters in the solution. We also tried to 
find exterior solutions corresponding to the Kerr solution in general relativity and have not 
yet succeeded in it. Note that the generic scalar-tensor- Weyl solution is not reduced to the 
Schwarzschild solution in the limit, tp — > 0. Therefore, we guess that the Kerr-like exterior 
solution must be reduced to the Tomimatsu-Sato-like solution in the limit, tp — > 0, where, by 
the term, the Kerr-like exterior solution, we mean the exterior solution such that it coincides 
with the Kerr solution when tp = and that it is reduced to the spherical exterior solution 
when the Kerr parameter (the angular momentum of the matter) vanishes. An example of the 
solution with the first property is obtained by using the metric function, ip, corresponding to 
the Kerr solution. However, the resulting solution does not satisfy the second property. We 
think that generating methods of the exact solutions (e.g. the Newman- Janis algorithm Jl3|) 
may play an important role in obtaining further exterior solutions. 

It is important to note that one have to specify an explicit form of the coupling function, 
A(tp), in order to obtain the complete information, especially the global geometry, of the 
exterior solutions in the scalar-tensor theories. We have not given any discussion on A(ip), 
and, instead, we have examined local geometrical properties of the solutions by taking account 
of the conformally invariant nature of null geodesies. We have evaluated the Ricci and Weyl 
deriving terms, TZ and F, in the scalar-tensor- Weyl solution and have found that local behaviors 
of null geodesies significantly depend on whether the parameter, 8, which corresponds to b/a 
in the spherical exterior solution with Q = 0, is greater than 1/2 or not. When S > 1/2, the 
behaviors are similar to those in the Schwarzschild spacetime. That is, the circular orbit exists, 
and the Weyl driving term, F, is strictly positive. When 8 < 1/2, on the contrary, there is no 



13 



circular orbit, and F changes its sign at xf- In the later case, one may expect to observe a 
gravitational lensing anomaly for images of distant sources appearing near \f- A ratio, 

C(h)= lim C+ ~ C ~ , (6.1) 

is a measure of the image deformation. In the Schwarzschild solution, C is a positive, decreasing 
function of h. However, in the spherical exterior solution with S < 1/2, C may become a 
negative increasing function of h near xf- Further quantitative studies of the solutions, C + 
and C_, and their astronomical applications, especially the implication to gravitational lens 
effects including their observational detectability, will be discussed in the forthcoming paper. 



A Summary of the scalar-tensor theories 

We shall consider the simplest scalar-tensor theory [2-5,14]. In this theory, gravitational 
interactions are mediated by a tensor field, g livi and a scalar field, 0. Hereafter, a symbol, 
denotes quantities or derivatives associated with g^ v . An action of the theory is the following: 



S- — 



4>R - ^r-9^j>J,v 



v /^rf 4 a; + 5 mattor [l' m ,.g^], (A.l) 



where u>(4>) is a dimcnsionlcss arbitrary function of 0, \fr m represents matter fields, and Setter 
is an action of the matter fields. The scalar field, 0, plays a role of an effective gravitational 
constant as G <~ 1/0. Varying the action by the tensor field, g^ v , and the scalar field, 0, yields, 
respectively, the following field equations: 

8tt^ , w(0) ( x 2 l ft Si , 2 \ , 1 



= -j- Tp„ + $ ^ v - -g^sTP^^p + T (V„V v - g„ v n4>), (A.2) 

6i - (A - 3) 

Now we perform the following conformal transformation to the unphysical frame (the Ein- 
stein frame) with the metric, g: 

= ( A - 4 ) 

such that 

GA 2 ( V ) = \, (A.5) 

where G is a bare gravitational constant, and A(ip) is referred to as a coupling function. Then 
the action is rewritten as 

S = ^ J(R- 2g^^.^. v )^—gd A x + S mattor [* m , A 2 (<p)g„ v ], (A.6) 
where the scalar field, ip, is defined by 

\ dtp J 3 + 2^(0) 
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Varying the action by g M „ and <p yields, respectively, 



1 



6 M „ = SirGT^ + 2 I cp^cp^ - ^g,i V g al3 <p, a (p,pj , 
Utp = -4TrGa{ip)T, 
where T^ v represents a energy-momentum tensor with respect to <?p„ defined by 
T ^ _ 2_8S mMtcr [*„,A 2 (<p) gilv ] = a6 ^ 



V~9 dg^u 
The conservation law of T^ v is given by 

^ V T» = a{<p)TV >. 

For the Maxwell field, one has 

TP U — — IF F a — -a F nF afi 



(A.8) 
(A.9) 

(A.10) 
(A.ll) 
(A.12) 



Further properties of the Maxwell field under the conformal transformation are found in Ref . | 
In this paper, we adopt the unit, G = 1. 



B Analytic solutions to the optical scalar equations 
B.l Static, spherically symmetric spacetime 



In this Appendix, we summarize the analytic results obtained by Dyer |l2|| . A metric of 
the static, spherically symmetric spacetime is given as 



ds 2 = ~e 2C dt 2 + e 2A dr 2 + e 2B dtf, 



(B.l) 



where A, B and C are functions of r. The null tangent vector, to a null geodesic in this 
spacetime is obtained as 



k° = c~ 2C , k 2 = 0, k' 6 = he 



-2D 



k 1 = ±e- ( - A+c Wl-h 2 c 2 (c- B \ 



(B.2) 



where a constant, h, is an impact parameter, and we assume that the geodesic is on the 9 = ir/2 
plane without loss of generality. The complex null vector, is obtained as 



t° = 



t L = 



-c 



2^2 
e~ A r 
2V2 



rS+S-e B + H{Sl + S 2 _) 



~(Sl + S 2 _)c B + 2HS+S- 
h 



t 2 = -^e" J 

, Hh 
V2 



-2B 



(B.3) 



■c B - c dr, 
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where a prime denotes a derivative with respect to r. Then the Ricci and Weyl driving terms, 
1Z and F, in (5.6) are evaluated as 

U + F = C -^ A + C ) (B" + {B'f - B'C - B'A') 

- h 2 e- 2( - A+B ) (C" + (C'f - C'B 1 - C'A') , 

Tl-F = c - 2 ( A +°) (B" + {B'f - B'C - B'A') (R4) 

- h 2 e^ A+B ^ (c 2< - A - B ^ + B" - A'B') . 

Since F is real, one can let a be real without loss of generality. Dyer has introduced the new 
optical scalars, C±, defined by 



— lnC± =9±a, 
d\ 



and has obtained the following equations: 

d 2 C ± 



d\ 2 



(TZ±F)C±. 



Dyer has obtained general solutions to (B.6) as 



C_ = C°_e B sin [ h 



e 

A+C 



2 LI 



P A+B+C | 

^dr +D° + \, 



h 2 e 2C }- 



. " =dr + D° . . 

^ C 2B Z. h 2 e 2C 1 



where C± and D± are integration constants. 



(B.5) 



(B.6) 



(B.7) 



B.2 Static, axisymmetric spacetime 



A metric of the static, axisymmetric spacetime is given by 

ds 2 = -e 2a dt 2 + e 2f3 dx 2 + e^dy 2 + e 2fi d(f> 2 , 



(B. 



where metric functions, a, (3, 7 and [i, are functions of x and y. The tangent vector, to a 
null geodesic on the y = plate is obtained as 



-2a 



, k 2 = 0, fc 3 = he- 2tl , 



k 1 = ±e~ l3 Ve- 2a - h 2 e- 2 », 
where a constant, h, is an impact parameter. The complex null vector, t^, is obtained as 



(B.9) 



1 



2V2 
1 

2V2 



-e> i S + S-+H(Sl + S 2 _) 



-(Sl + S 2 _)e t * + 2HS+S_ 
h 



c-«, t 2 = ^c- 
1 



t 3 = -^Hhe~ 211 , 



V2 



(B.10) 



S± = Ve~ a ± hc-», H=-- 

h 



S+S- dx 



dx. 
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Figure 2: The Weyl and Ricci deriving terms, F and 1Z, are shown, where (A, S, h) = 
(1.5,0.3,2,0) in (a) and (1.5,0.7,2.0) in (b). The vertical and horizontal axes denote, re- 
spectively, the deriving terms and x. 
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Scalar-Tensor- Weyl ( A , 6 ) 



A — 1 
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A ► 8 



5 — ► 1 



Weyl (5) 



5^ 1 



Schwarzschild 



Figure 1: A schematic picture of series of the scalar-tensor- Weyl solutions. 
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Figure 2: The Weyl and Ricci deriving terms, F and 1Z, are shown, where (A, <5, h) 
(1.5,0.3,2,0) in (a) and (1.5,0.7,2.0) in (6). The vertical and horizontal axes denote, 
spectively, the deriving terms and x. 



